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The Cygnus region hosts one of the most remarkable star-forming regions in the Milky Way. In¬ 
deed, the total mass in molecular gas of the Cygnus X complex exceeds 10 times the total mass 
of all other nearby star-forming regions. Surveys at all wavelengths, from radio to y-rays, re¬ 
veal that Cygnus contains such a wealth and variety of sources—supernova remnants (SNRs), 
pulsars, pulsar wind nebulae (PWNe), Hii regions, Wolf-Rayet binaries, OB associations, micro¬ 
quasars, dense molecular clouds and superbubbles—as to practically be a galaxy in microcosm. 
The y-ray observations along reveal a wealth of intriguing sources at energies between 1 GeV 
and tens of TeV. However, a complete understanding of the physical phenomena producing this 
y-ray emission first requires us to disentangle overlapping sources and reconcile discordant pic¬ 
tures at different energies. This task is made more challenging by the limited angular resolution 
of instruments such as the Fermi Large Area Telescope, ARGO-YBJ, and HAWC and the limited 
sensitivity and held of view of current imaging atmospheric Cherenkov telescopes (lACTs). The 
Cherenkov Telescope Array (CTA), with its improved angular resolution, large held of view, and 
order of magnitude gain in sensitivity over current lACTs, has the potential to hnally create a co¬ 
herent and well-resolved picture of the Cygnus region between a few tens of GeV and a hundred 
TeV. We describe a proposed strategy to study the Cygnus region using CTA data, which combines 
a survey of the whole region at 65° < I < 85° and —3.5° < b < 3.5° with deeper observations of 
two sub-regions that host rich groups of known y-ray sources. 
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1. Introduction 

Among the star-forming regions of the Milky Way, the Cygnus region stands out as a superb 
laboratory for cosmic ray production, acceleration and transport. Surveys of the region at all wave¬ 
lengths, from radio to 7 -rays, reveal a galaxy in microcosm, filled with supernova remnants (SNRs), 
pulsars, pulsar wind nebulae (PWNe), H n regions, Wolf-Rayet binaries, many OB associations, mi¬ 
croquasars, dense molecular clouds and superbubbles, and even a 7 -ray source that testifies to the 
presence of cosmic ray sources within superbubbles [1]. This last lies in the nearby (J 1.5 kpc) 
Cygnus X complex, ten times as rich in molecular gas as all the other nearby star-forming regions 
[2, 3]. Multi-messenger data from IceCube and AMANDA data hints at a neutrino signal from 
the Cygnus region, which, if confirmed, would be the unambiguous signature of powerful proton 
acceleration above tens of TeV [4]. 

Cygnus also boasts diffuse VHE 7 -ray emission in excess of that predicted by modeling the 
interactions of the large-scale cosmic-ray population with interstellar matter and radiation fields 
in the Galaxy [1, 5]. This diffuse excess may indicate either a local overabundance of freshly- 
injected cosmic rays or an unresolved population of faint very high energy (VHE; E > 100 GeV) 
7 -ray-bright sources. Einally, the Cygnus region hosts two of the most remarkable high-mass X- 
ray binaries known: Cyg X-3 and Cyg X-1. Very sensitive observations of the Cygnus binaries in 
7 rays could permit us to detect and probe the origin of VHE emission in microquasars and thus 
access the physics of the jets and their interaction with the environment. [ 6 , 7]. 

1.1 Current Status 

In Cygnus, the brightest portion of the Northern 7 -ray sky, we find a wealth of intriguing 7 -ray 
sources at energies between 1 GeV and tens of TeV. However, the picture provided by the various 7 - 
ray observatories— 7 -ray satellites such the Eermi Earge Area Telescope {Fermi EAT ) and AGIEE 
at GeV energies, imaging atmospheric Cherenkov telescopes (lACTs) such as HEGRA, MAGIC 
and VERITAS at hundreds of GeV, and Milagro and ARGO-YBJ [ 8 ] at TeV energies—is complex 
and at times discordant. Beginning at 1 GeV Fermi EAT detects an extended excess of hard 7 -ray 
emission, the so-called Eermi cocoon [1]. ARGO-YBJ and Milagro have detected sources that they 
associate with the cocoon at higher energies [9]. AGIEE also reports a tentative detection of the 
X-ray persistent stellar mass black-hole binary Cyg X-1 in a flaring state [10]. 

Surveys and targeted snapshots by lACTs present a different view. In addition to MAGIC’s 
reported detection of Cyg X-1 during a flare, they also identify a set of more compact extended 
sources: TeV J2032-I-4130, the first unidentified and extended source ever discovered at very high 
energies, VER J2019-I-407, likely linked to radio-. X-ray-, and 7 -ray-bright SNR G78.2 -|-2.1, and 
VER J2019-I-37, linked to MGRO J2019-I-37 and one of the hardest-spectrum sources ever detected 
[11, 12, 13, 14]. Only upper limits could be set for the other well-known microquasar in the region, 
Cyg X-3, which is located in a much denser environment that Cyg X-1 [15, 16]. Conspicuous by 
its absence is any lACT detection of the Cygnus cocoon. 

We must disentangle these overlapping 7 -ray sources if we are to definitely pinpoint the phys¬ 
ical phenomena behind them. This task is made more challenging by the limited angular resolution 
of instruments such as the Fermi EAT , ARGO-YBJ, and HAWC and the limited sensitivity and 
field of view of both past and current lACTs. Establishing whether freshly-accelerated cosmic rays 
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or unresolved faint sources [17, 18] produce the bright VHE 7 -ray diffuse emission detected by 
MILAGRO [19, 20] also requires more sensitive instruments with finer angular resolution. 

2. Role of CTA 

CTA has the capabilities needed for a fruitful study of Cygnus [21]. A factor of ten improved 
sensitivity over current lACTs provides the statistics needed for detailed morphological and spec¬ 
tral studies as well as studies of transient phenomena from high-mass X-ray binaries. The large 
(7 — 8 °) field of view (FoV) is well-suited to surveys and boosts sensitivity to extended and diffuse 
sources. Finally, CTA’s angular resolution far surpasses that of current 7 -ray observatories, includ¬ 
ing HAWC, between a few hundred GeV and ~ 100 TeV, which will help separate out overlapping 
sources and pinpoint low-energy counterparts. CTA’s planned program of Key Science Programs, 
or KSPs, provides for a survey of the region 65° <l < 85° and —3.5° <b <3.5° down to 2.7 mCrab 
under the aegis of the Galactic Plane Survey KSP The Star-forming Regions KSP will augment this 
survey with deep observations of two sub-regions that host the richest and most interesting groups 
of known 7 -ray sources: the region surrounding Cygnus OB2/Fermi cocoon and that surrounding 
Cygnus OBl/MGRO J2019-I-37. We discuss the projected impact of these observations below. 

2.1 Cygnus OB2 and the Fermi Cocoon 

Source confusion presents a major challenge when trying to pinpoint the source of the ener¬ 
getic particles flooding Cygnus-X. The favored option is the collective action of star-wind shocks 
in the stellar clusters, but contributions from SNR G78.2-I-2.1 or a PWNe cannot be completely 
ruled out. Fikewise, the energy spectra of the Fermi FAT cocoon, ARGO J2031-I-4157, and 
MGRO J2031-1-41 (hard 7 - ray excesses of roughly ~ 2° extension co-located between Cyg OB2 
and SNR G78.2-1-2.1) connect relatively smoothly, bolstering the interpretation that these are detec¬ 
tions of the same source at three different energies. Yet the coarser angular resolution of ARGO- 
YBJ and Milagro means that these instruments must rely on lACT measurements to correct the 
VHE spectrum for contributions both VER J2019-1-407 and TeV J2032-I-4130, which lies near Cyg 
OB2 and surrounds a 7 - ray pulsar that may itself be part of a binary system [22]. This incurs 
an inflated systematic uncertainty due to cross-calibration issues [9] and is further exacerbated by 
the non-detection of the cocoon by current lACTs, likely to due background-subtraction effects 
connected to the small field of view. A more self-consistent treatment of the problem requires an 
instrument, like CTA, that can both detect and resolve all sources in the region. 

CTA may also shed light on the nature of VER J2019-I-407, which occupies a region of the 
northwest rim of G78.2-1-2.1 bright in radio and diffuse non-thermal X-rays [12] that is less than 
~ 16% of the area of the SNR. As figure 1 shows, the sources share a common power-law spectral 
index (E ~ 2.4) but differ widely in flux as well as angular extent. These two sources may reflect 
accelerated particle populations with different energy cutoffs in different parts of the remnant, with 
the bulk of the remnant having a much lower cutoff (e.g. the green curve in figure 1) or VER 
J2019-I-407 may be something else entirely, such as a PWN in the line of sight. CTA’s better 
angular resolution and large field of view will make it much easier to disentangle the contributions 
from the remnant and VER J2019-1-407 and assess both their mutual relationship and their effect on 
the measured spectrum of the cocoon. 
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Figure 1: Best-measured spectra above 10 GeV for all known sources in this region apart from the micro¬ 
quasar Cyg X-3 and the contributions from 7 -ray pulsars, whose low energy cutoffs make them insignificant 
at these energies. Dashed lines show plausible extrapolations of the fitted spectra. Two independent measure¬ 
ments by VERITAS (red points) and MAGIC (red solid line and bow tie) are shown for TeV J2032H-4130. 


This is home out by a dedicated simulation of a 100 hour observation of the Cyg OB2/cocoon 
region, processed through a binned likelihood analysis. These simulations were run at the Cyfronet 
computing facility of the University of Cracow using ctools for CTA and the instrument response 
function was based on a site at Tenerife. All known sources of VHE 7 -ray emission above 10 
GeV (the Fermi Cocoon, SNR G78.2-t2.1, VER J2019-t4007, TeV J2031-t4130) were included 
and all four are detected with high (> 18a) significance, with a number of detected 7 -ray counts 
which exceeds 5% of the residual CR background counts integrated over the source area. The 
count distribution as a function of energy (Eigure 2a) can be reliably used to determine arcmin- 
level positions and extensions of the sources with few % statistical uncertainties. The spectra of the 
sources are also reliably reconstructed at the few % statistical level, including all cutoff energies 
below 10 TeV (Eigure 2b). Precisely constraining cutoff energies above 10 TeV will require deeper 
exposures. 

This work may benefit in future from a more advanced likelihood technique, which uses a 
7 /hadron classification parameter to gain additional separation power between photons and cosmic 
rays. Such a technique improves overall sensitivity and compensates to some degree for influence 
of fhe cosmic ray background, which sfrongly impacfs sensifivify below 300 GeV and confribufes 
fo fhe fendency of fhe fil fo underesfimafe fhe extension of fhe more compacf sources. Nonefheless, 
we may be confidenl fhaf fhe CTA observafions of Cygnus will shed lighf on whefher fhe cocoon 
is a single source or fhe superposifion of multiple objecfs and clarify fhe relafionship befween 
VER 2019-1-407 and SNR G78.2-I-2.1. The excellenf specfral characterization will help esfablish 
whefher fhe 7 -ray emission is due primarily fo hadronic or lepfonic processes for fhese objecfs. 
The measured cutoff energies will play a critical role in esfablishing maximum particle acceleration 
energies in fhe former case and probing energy losses in exfreme environmenfs in fhe laffer. 
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Figure 2: Results from simulated deep-field observations of the Cygnus OKI!Fermi LAT cocoon region. 
Left; Simulated counts maps. Right; Reconstructed spectra of known sources. 

2.2 Cygnus OBI and MGRO J2019-I-37 

The Cygnus OBI region simulation considered known VHE y-ray emitters MGRO J2019-1-37, 
VER J2019-I-368, and VER J2016-I-371. In all cases the simulations extrapolated from the measured 
spatial and spectral properties of these sources. The observed spectrum of the MGRO 12019-1-37 
seems to be best fitted with power law with exponential cutoff model. The spectral index is found 
to be r = 2.0 and the cutoff energy Ec = TeV, with large error bars. No cutoff indication 

is given for VER J2016-I-371 and VER J2019-I-368. In order to understand whether CTA will be 
able to study the morphology of the Cyg OBI region and to resolve high energy cutoffs in the 
energy spectra of these sources, MGRO J2019-I-37, VER J2019-I-368 and VER J2016-I-371 were 
each simulated with cutoff energies equal to 10, 30 and 80 TeV. 

Spectra are extracted from the resulting simulated data using the same analysis as in Sec. 
2.1; the spectra are shown in figure 3. The simulafed cutoff energies for MGRO J2019-I-37 are 
reasonably well-reproduced when fil by power law model wifh an exponenfial cutoff; 13.2 ±0.6 
TeV, 29.3 ±1.9 TeV and 77.5 ±9.5 TeV. We anficipafe fhaf CTA has greaf pofenfial to nail down 
fhe specfral culoff of fhis source af high energies, where fhe VHE y-ray emission is almost free 
of cosmic ray background. The spectra of VER 12019-1-368 and VER 12016-1-371, however, fit 
better to a power law model. Notably, the spectral analysis using an exponential cutoff failed for 
VER 12019-1-368 when the assumed cutoff energy was chosen at 10 TeV, while the fitting errors 
exceeded 100% for the 30 and 80 TeV cases. 

2.3 Physics of binary systems 

Cyg X-3, identified as fhe firsf GeV y-ray-emiffing microquasar, has opened up a new window 
on fhe non-fhermal processes associafed wifh relafivisfic jefs from accreting black holes. The Fermi 
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Figure 3: Spectra obtained from fits to simulated data for MGRO J2019+37 (top left), VER J2019+368 (top 
right) and VER J2016+371 (bottom) in the Cygnus OBI region, each with three different assumed energy 
cutoffs (red, 10 TeV; green, 30 TeV; blue, 80 TeV). 


LAT orbital light curve of Cyg X-3 suggests that GeV emission may originate from a region well- 
separated from the binary system. It could be either due to the jets (via internal shocks) or from 
a re-collimation shock when the jets interact with the dense wind from the Wolf-Rayet stellar 
companion. At this stage, the y-ray spectrum of Cyg X-3 is still unconstrained at VHE. 

Extrapolating the spectrum of the Cyg X-3 flare seen by Fermi EAT [23] at energies > 30 
GeV we estimate that a significant detection will take 200 min of observation time with 10 hours of 
exposure needed for the determination of the full spectral shape. It is therefore plausible that CTA’s 
planned observations of Cygnus could trap a flare from either Cyg X-1 or Cygnus X-3. Persistent 
emission is reported with Eermi in the case of Cyg X-3 [24] with a flux of ~ 0.5x10 ^phcm ^s ^ 
in the 0.1-10 GeV range. Without further spectral information extrapolation into the energy range 
probed by CTA is very difficult. However, the few mCrab sensitivity promised by the CTA KSPs at 
the position of Cyg X-3 should strongly constrain both the level and nature of the high-energy y-ray 
emission in the event of a non-detection. CTA should be able to constrain the maximum energy 
reachable by this cosmic accelerator [25] and potentially localize the origin of the VHE emission. 

Cyg X-1: INTEGRAE’s recent detection [26] of strongly polarized hard X-ray emission from 
Cyg X-1 supports the existence of high-energy particles in this source. A CTA detection of Cyg 
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X-1 would help to bridge the gap with stellar mass black holes in low mass X-ray binaries, which 
undergo regular outburst with a specific phase of transient jet emission. If we assume a Cyg X- 
1 flare with similar properties to that reported MAGIC {dN/{dAdtdE) = /ITeV)^^'^) 

flare [13], we find from back-of-the-envelope estimates that a 5-a detection above 30 GeV can 
be achieved in 2 minutes. 10 minutes will be required to measure the spectral shape. Persistent 
emission from Cyg X-1 has been detected at 1 MeV using INTEGRAL [26] but no corresponding 
detection has been reported by Eenni EAT. A naive extrapolation of the INTEGRAL detection to 
CTA energies predicts emission at the mCrab level, roughly comparable to the sensitivity of the 
planned CTA GPS at this position. 

2.4 The Milagro Diffuse Emission 

We investigated the number intensity relation for HESS-detected VHE y-ray sources above 20 
mCrab in the region of longitude 250° <l < 65° and latitude —3.5° <b < 3.5° [27]. Extrapolating 
the number intensity curve down to 3-4 mCrab (the sensitivity planned for the Northern GPS) 
about 600-500 sources could be detected by CTA from the HESS region. Rescaling this number 
for Cygnus as in [17], faint sources resolved by CTA could account for about 40 per cent of the 
diffuse emission measured by Milagro. CTA can therefore better separate the true diffuse TeV 
y-ray emission from unresolved sources and shed light on its nature. 
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